Single-atom lasing induced atomic self trapping 
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We study motion and field dynamics of a single-atom laser consisting of a single incoherently 
pumped free atom moving in an optical high- Q resonator. For sufficient pumping, the system starts 
lasing whenever the atom is close to a field antinode. If the field mode eigenfrequency is larger than 
the atomic transition frequency, the generated laser light attracts the atom to the field antinode 
and cools its motion. Using quantum Monte Carlo wave function simulations, we investigate this 
coupled atom-field dynamics including photon recoil and cavity decay. In the regime of strong 
coupling, the generated field shows strong nonclassical features like photon antibunching, and the 
atom is spatially confined and cooled to sub-Doppler temperatures. 

PACS numbers: 32.80.Pj, 42.50.Vk, 42.50.Lc 



Optical cavity QED experiments using high finesse res- 
onators and ultracold neutral atoms have seen tremen- 
dous progress towards larger coupling strength and in- 
teraction time in the past decade, becoming a fruitful 
test ground of quantum theory 0, M EI S HI O ne par- 
ticular goal of cavity QED is a single-atom, single-mode 
laser. In analogy to the micromaser, a one atom laser 
using an atomic beam crossing an optical resonator has 
been realized some time ago [6j- However, although effec- 
tively only a single or a few atoms are present at a time, 
the short interaction times make this still a many particle 
device. Very recently important progress towards single 
particle lasing has been achieved by placing a small ion 
trap into a high- Q cavity 0, ■ This allows very long 
interaction times, but technical limits set by the ion trap 
still prevent to reach threshold. 

In all optical setups the atoms are extremely cold at the 
beginning, but the light forces induced by the cavity field 
and pumping itself lead to fast heating. Even less than a 
single photon on average induces significant dissipation 
Q and limits interaction times. By a proper choice of pa- 
rameters, one can minimize this heating caused by light 
forces || or even use cavity induced cooling forces 0, [H3 
to enlarge the trapping times. Alternatively, the addition 
of an extra dipole trapping potential led to much larger 
interaction times and enabled the first realization 
of single atom lasing in the strong coupling regime 
Here, steady state light output generated from one atom 
for almost a second was achieved from a three photon 
Raman type gain scheme. 

In this Letter we go conceptually one step beyond and 
assume that no coherent field is applied to the atom or 
the mode and no extra atom trap is present. The cavity 
field is entirely generated by gain from the incoherently 
pumped atom. Neglecting atomic motion, such single 
atom laser models have been widely used in fundamental 
studies of quantum laser theory [12(. They constitute 
light sources with unique properties. For example, one 
is able to sustain a stationary single photon Fock state 
or generate highly sub-Poissonian output 13) ■ Note that 
the same single atom is present during the whole time 



here. 

Extending these models, we include the atomic mo- 
tion governed by the light forces of the laser field created 
by the atom itself. As expected, the proper accounting 
of the cavity field forces implies significant modifications 
of the system dynamics In general, the atomic mo- 
tional variables and the dynamics of the internal vari- 
ables are coupled and correlated or even entangled. The 
atom moves under the influence of the dipole force which 
depends on the mode intensity. This intensity in turn de- 
pends on the atomic position in a highly nonlinear way 
as the atom itself is the gain medium. 

From a first guess, one could already expect some self- 
trapping effect if one chooses the parameters such that 
lasing only starts when the atom is close to a field antin- 
ode. The laser field then could generate an attractive 
potential keeping the atom in the vicinity of this antin- 
ode. Naively, this requires the lasing mode to be red 
detuned relative to the atomic transition frequency. In 
this case, though, the extra energy from the atomic tran- 
sition is likely to be converted into kinetic energy heating 
the atom. Similarly, for a blue detuned cavity mode, one 
expects cooling but repulsion from the antinodes. How- 
ever, the situation is a bit more complicated. As gain 
requires atomic inversion, an atom in steady state can 
still be a high field seeker for blue detuning, and the 
laser frequency is itself a dynamical quantity. Thus, we 
can still hope to find parameters where trapping, cool- 
ing, and lasing coincide. Of course, heating through 
spontaneous emission and dipole fluctuations will also be 
present. This makes the total dynamics hard to guess, 
which motivated us to study the problem in more detail. 

Let us now define our model as simple as possible still 
containing the essential physics. For this we restrict our- 
selves to a two-level atom moving in a single strongly 
coupled cavity mode in one dimension. As in well proven 
approaches developed in the early days of laser theory 
[12| , incoherent pumping can be consistently modelled by 
inverse spontaneous emission at rate 26. Following stan- 
dard procedures of quantum optics, we can derive the 
master equation p = Cp governing the time evolution of 
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field and atom including atomic spontaneous emission at In a frame rotating at the mode frequency lo c we get: 
rate 27 and cavity decay at rate 2k. 
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Here \g) (|e)) and a denote the atomic ground (excited) 
state and the field annihilation operator, respectively. 
The cavity mode with mode function cos(fca;) and fre- 
quency lu c is detuned from the atomic transition fre- 
quency L0 a by A = lu c - uj a . 

In a first step we look at the steady state of the sys- 
tem for fixed atomic position x which enters the equations 
only via the coupling strength g cos(fcx) . In Fig.^we plot 
the photon number n (solid line), its scaled uncertainty 
A n /n (dashed), and the atomic upper state population 
(dash-dotted) as a function of g for fixed pump strength 
5 = 75k = 7.57 for large atom-field detuning A = 250k. 
As expected, the photon number depends upon g in a 
nonlinear way and the system starts lasing only for suffi- 
ciently large g when the atom is close to a field antinodc. 
With growing photon number (n > 1) the correspond- 
ing spectrum shown in the little inserts of Fig. ^ is blue 
shifted from the atomic resonance and acquires a width 
below the empty cavity linewidth. Hence, we can hope 
for stable trapping close to the antinode combined with 
lasing at these parameters. 

Let us include atomic motion and check. Since the 
number of degrees of freedom is rather large now, we 
use a Monte Carlo wave function simulation technique 
to numerically approximate the solution of the master 




FIG. 1: Average photon number n, scaled uncertainty An/n, 
and upper state population as function of gjn. The param- 
eters are (7,(5, A) = (10k, 75k, 250k). The little inserts show 
the emitted light spectrum for g = (5, 25, 45)k. 
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FIG. 2: Average photon number n and Mandel Q parameter 
for the lasing startup phase. The parameters are (j,5,g) = 
(10k, 60k, 50k) where the cavity field is detuned from the 
atomic transition by A = 250k. 



equation Averaging gives an approximate density 

operator p(t) while analyzing individual trajectories pro- 
vides extra insight in the microscopical dynamics. 

We briefly review the method now. Defining a super- 
operator Sp = J2i CipC[ based on collapse operators C; 
related to the quantum jumps, we get single trajectories 
by propagation with a non-hcrmitian Hamiltonian Heff 



(£ — S)p = —iHcftp + ipHeW, 



(2) 



in between stochastically occurring jumps. The jump 
probabilities within the interval [t, t + At] are given by 
Pi(t) = (ip{t)\ClCi\ip(t))At. One first propagates \tp(t)) 
using HcS for one time step At 



\ip(t + At)) = exp(-iHesAt)\ip(t)) 



(3) 



an calculates p,. Using random numbers <G [0, 1] and 
comparing them to Pi(t), one then decides on the occur- 
rence of collapses — > Ci\ip) at the end of At. The 
jumps here include atomic decay (^/2r/\g) (e\) , cavity de- 
cay (V2ko), and pump events (V2S\e)(g\). 

Explicitly the effective Hamiltonian Hefi reads (h = 1) 

HoS — — iS\g)(g\ — (A + Z7)|e)(e| — ina^a 

- igcos(kx) (\e)(g\a - a^\g)(e\) . (4) 
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FIG. 3: Mean square velocity in units of the Doppler tem- 
perature Td for A = 250k;. The other parameters are 
(7, 8, g) = (10k, 60k, 50k), Td is indicated by the dashed 
line. For negative detuning the atom gets heated up, see inset 
where A = —250k. 

As the atomic temperature is well above the recoil 
limit, we treat the external atomic variables x and p clas- 
sically [T^|. Thus we have A = p/m and p = F to add to 
our equations. The force F is given by 

dUcff 

F = (i>\ K— W) = ~ il% k sin kx(ip\e)(g\a\tp)+h.c. (5) 

ox 

Finally, wc have to truncate the Hilbert space H = 
span{|<7), |e}} ® tip (TCf is the Fock-space of the mode) 
at photon number N. Our state space thus contains the 
ground state \g, 0) and N manifolds separated by u> a with 
states \g, n) and |e, n — 1), n = 1...N. Any wave function 
is determined by 2N + 1 complex numbers, 

JV 

= .9o|.9,0) + ^2 (9n\g,n) +e n \e,n- 1)) , (6) 

n=l 

where only 2N are independent due to normalization. 
Note that in contrast to most previous treatments of light 
forces in a cavity, we cannot adiabatically eliminate the 
excited state since we have to deal with an inverted atom. 

Let us now examine stochastic simulation averages for 
typical cases. After turning on the pump, the mean pho- 
ton number n increases rapidly and reaches a constant 
value within a few cavity relaxation times (Fig. El . At the 
same time the atom gets trapped and remains captured 
in the potential of the light it has generated. In addtition 
we show the Mandel Q parameter (Q = (n 2 — rl 2 )/ (n) — 1) 
as a measure of the field intensity noise. 

Fig. |21 shows the atomic kinetic energy averaged over 
4000 realizations in units of the Doppler temperature 
Td = ^7/2 for the parameters of Fig. [3 where the mean 
photon number reaches n = 3.8. Starting at a very low 
velocity, the atom gets cooled (heated) for a wide range 
of positive (negative) detunings A. The inset shows the 
analogous result for A = —250k. This agrees with a sim- 
ple energy conservation argument. The atom gains the 



FIG. 4: Position distribution of an atom trapped at an antin- 
ode for the parameters of Fig. [3] (solid line) and for a larger 
pumping rate S = 80k (dashed line). 

internal energy Tiuj a from a pump event. Subsequently 
it loses the energy Tiu) c by stimulated emission into the 
cavity mode. The energy mismatch is compensated by 
the atomic center of mass motion which gets damped if 
lu c — uj a > 0. This argument is confirmed by a closer 
look at basic absorption and emission processes includ- 
ing Doppler shift from the atomic motion. 

Naively one would guess that the atom is expelled from 
the interaction region of a blue shifted light field. How- 
ever, this is not true for a partially inverted atom which 
still can be trapped near the mode antinodes. This fea- 
ture is essential to allow steady state operation of our 
single-atom laser. The solid line in Fig. 2] shows the po- 
sition distribution of an atom trapped after the cooling 
process of Fig. [3] Clearly, the atom spends most of the 
time in the vicinity of an antinode with a mean square 
distance of x 2 = (0.1A) 2 . This value decreases further 
for larger photon numbers. For n f=a 12 (dashed line) we 
get x 2 = (0.07A) 2 for a pump strength of S = 80k. The 
non- vanishing probability near x « ±A/4 shows some re- 
maining hopping of the atom between different trapping 
sites. 

So far we calculated ensemble averages to investigate 
the motional characteristics. It is now quite interest- 
ing to look at individual trajectories \ip(t)) to visualize 
the microscopic dynamics. Fig. OJt depicts a typical ex- 
ample of the photon number time evolution (n) for a 
trapped atom for ten cavity relaxation times. Starting 
in a certain state |e,n), the system gets entangled with 
\g, n + 1) owning to the atom-cavity coupling. The inco- 
herent pumping projects it into the state |e, n+ 1) and 
creates a Fock state with n + 1 photons. Unless a pho- 
ton leaks out of the cavity or another one is created, the 
photonic state remains nearly unchanged for the follow- 
ing reasons. First, we chose a large detuning A = 250k in 
order to keep the atom trapped at antinodes as good as 
possible. Second, the atom undergoes several jumps due 
to alternating spontaneous emission and pump events. 
After each of these cycles the system is again projected 
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FIG. 5: a Photon number expectation value (n) for a sin- 
gle trajectory of a trapped atom within ten cavity relaxation 
times. The parameters are as in Fig. 01 The dashed line 
indicates the mean photon number, b Corresponding single- 
trajectory q factor. 
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FIG. 6: Mean photon number and Mandel parameter vs 
pumping strength S. 



into the initial state \e,n + 1), which gives the cascade 
shape of the photon number time evolution. Notice that 
the pumping process creates no coherence, as it would 
be e.g. in a lambda type system such as in [T^]. Here ac 
Stark splitting of dressed states does not effect the prob- 
ability of exciting the system which is solely given by the 
atomic ground state population and, hence, the system 
can be excited into states with higher photon numbers 
more easily. Nevertheless, for each trajectory the cavity 
field remains close to a Fock state with photon numbers 
varying around the mean value indicated by the dashed 
line. This behavior is demonstrated in Fig. \Sjp where 
the single trajectory Mandel factor - which describes the 
actual state of the cavity field - attains values close to 
q = — 1 as long as one photon is present at least. This 
is due to the fact that in each trajectory we know the 



number of pump and decay events. Statistical averaging 
over the pump events washes out this feature. 

Let us now come back to ensemble averages of field 
properties. Fig. depicts n and Q function of 5 
with (A, g) = (200k, 100k) for 7 = 0. In general the field 
intensity is strongly fluctuating. These fluctuations are 
even more pronounced than for an atom at rest and in 
turn increase motional heating. Only for larger photon 
numbers the Mandel Q parameter drops down to zero as 
for a coherent state (lasing). 

In summary light forces significantly influence the dy- 
namics of a single-atom laser. Surprisingly, for blue 
atom-field detuning several effects work together in a 
favorable way to facilitate steady state lasing in con- 
junction with mechanical cooling and trapping. Steady 
state temperatures below the Doppler limit are possible 
in spite of enhanced momentum diffusion due to fluc- 
tuations of the photon number. The light field shows 
nonclassical features and approximates a coherent state 
with a coherence time beyond the cavity lifetime only far 
above threshold. Although the description of the atom is 
oversimplified here, we believe that our central findings 
like enhanced trapping and self cooling with lasing still 
should be present in a real system. 

This work was supported by the Austrian FWF un- 
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